Nonlinear cotunneling through an artificial molecule 
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We study electron transport through a system of two lateral quantum dots coupled in series. We 
consider the case of weak coupling to the leads and a bias point in the Coulomb blockade. After 
a generalized Schrieffer- Wolff transformation, cotunneling through this system is described using 
methods from lowest-order perturbation theory. We study the system for arbitrary bias voltages 
below the Coulomb energy. We observe a rich, non-monotonic behaviour of the stationary current 
depending on the internal degrees of freedom. In particular, it turns out that at fixed transport 
voltage, the current through the system is largest at weak-to-intermediate inter-dot coupling. 

PACS numbers: 73.63.Kv, 73.23.Hk, 72.10.-d, 03.67.Lx 
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Quantum dots are prototype systems for studying 
the properties of discrete levels embedded in a solid- 
state environment. Single dots ("artificial atoms"—) 
can be coupled through quantum point contacts, lead- 
ing to "artificial molecules". Indeed it has been shown 
experimentallySi^, that the eigenstates of double dot 
systems are coherent molecular superpositions of single 
dot (atomic) states. Unlike real molecules, these dots 
are readily contacted and tunable in situ, making them 
a natural test-bed for molecular transport. Double dots 
have also been proposed as charge quantum bits^. 

This raises the question, which information on the en- 
ergy spectrum and the wavefunctions of the dot can be 
probed by transport measurements. This is only pos- 
sible if artifacts induced by the coupling to the leads 
can be sorted out and the double-dot is disturbed as 
little as possible. This is the case, when the coupling 
to the outside leads is weak (see Fig. and the gates 
are tuned to the Coulomb blockade regime^i^. In that 
regime, only states with a fixed number of electrons are 
energetically permissible and hence sequential tunneling 
is suppressed. The leading transport mechanism in this 
case is cotunneling^^, the coherent transfer of two elec- 
trons via virtual levels in the dots. Our work stands be- 
tween studies focusing on sequential tunneling^ and work 
on linear response in the Kondo regimeiS. The proper- 
ties of cotunneling currents as a spectroscopic tool for 
the spectrum of quantum dot system have recently been 
studied in exquisitely controlled experiments on systems 
similar to ours^ti^. 

In this paper, we analyze a serial configuration of lat- 
eral quantum dots in the cotunneling regime. We study 
finite voltages up to the order of the charging energy i.e. 
do not restrict ourselves to linear response. We find a rich 
nonmonotonic structure in the current as a function of 
the dot parameters. In particular, we find a pronounced 
crossover indicating the opening of an inelastic transport 
channel, which leads to the surprising result, that a too 
strong interdot coupling actually inhibits charge trans- 
port. We analyze the infiuence of the asymmetry of the 
dots on the current. 

In the Coulomb blockade regime^^, the relevant 
Hilbert space is spanned by four basis states i,j G 



{0, 1}, which denotes i and j additional electrons (as 
compared to an appropriate neutral state) on the left 
and right dot respectively. We study the situation where 
the gate voltages of the single dots are very close to each 
other and the inter-dot coupling is, although appreciable, 
much smaller than the single dot addition energy. Thus, 
the subspace spanned by the two states |1,0) and |0, 1) 
is energetically most favorable. The next closest states 
\vo) = 1 0,0) and \v2) = |1,1) are outside the transport 
window and serve as virtual statesii. States with higher 
dipolar moment are energetically even less favorable due 
to the high charging energy of the individual dots. 
The Hamiltonian of this system can be written as 
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Note, that the sum over dot states n only runs over the re- 
stricted Hilbert space described above. Hq describes the 
isolated double-dot {Hgys) and the leads {H^es), whereas 
the tunneling part Hi describes the coupling of each dot 
to its lead and will be treated as a perturbation, fii/r 
are the number operators counting additional electrons 
on either dot. The asymmetry energy Sas = (e; — er)/2 
describes half of the difference between the energy level 
for the additional electron in left dot (e;) and the cor- 
responding energy level in the right dot (Cr), which can 
be tuned through the gate voltages, e/j and are the 
charging energies towards the higher level \v2) and the 
lower level \vq) respectively. 7 is the tunable inter-dot 
coupling strength. The a'^^s and &*^'''^s denote electron 
creation/annihilation operators in the dots and leads. In 
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FIG. 1: Sketch of the considered artificial molecule, where 25 
is the level splitting and V the bias voltage. The coupling 
to the outside leads (hatched areas) is assumed to be small 
whereas the inter-dot coupling (dotted line) can be strong. 



Hi, the symbol tc represents the tunnel matrix element 
between the dots and the leads. It is independent of the 
energies in the double-dot system and the correspond- 
ing sequential tunneling rate fiV = 2TTtlN {e-p) should be 
small compared to the internal energies. N{£-p) is the 
density of states in the leads taken at the Fermi energy. 
We restrict our analysis to spin-polarized electrons, these 
can be polarized by an appropriate in-plane magnetic 
field. Fig. ^ shows a sketch of the system. In Fig. ^ 
V — /jR — iiL is the bias voltage between the two leads 
(hatched) and 26 = + 7^ is the level splitting in 

the molecular two-state system (TSS). 

Pursuing our aforementioned objective, we take the 
inter-dot coupling 7 into account to all orders by diag- 
onalizing -ffgys ^^nd transforming Hi into the new basis. 
Already now, there is no simple selection rule or symme- 
try of the coupling of the states to the leads any more. We 
want to use well-established tools of lowest order pertur- 
bation theory for both finding the density matrix of the 
system and evaluating the current. In order to capture 
cotunneling by this approach, we perform a Schrieffer- 
Wolff transformatiorki^ up to second order, i.e. we take 
into account all indirect transitions between arbitrary fi- 
nal and initial states of the dot which involve only a single 
intermediate state. This takes the transformed Hamilto- 
nian into the generic form 
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are Schrieffer- Wolff amplitudes and 
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where the H. '. 

k,k' 

c,d — ± denote the two molecular levels, a).'/'^ the as- 
sociated molecular operators and Y, Y' the position of 
the electrons involved in these processes. Due to the 
molecular nature of the double dot eigenstates, all the 
amplitudes are finite and composed out of a huge num- 
ber of contributions with no particular symmetry. The 
perturbation theory formula for this general case can be 
found e.g. in Ref. [3 Etnd is worked out in more detail 
in Refs.'w andlTsl In eq. we have taken matrix ele- 
ments in the double-dot eigenbasis only whereas we stick 
to second-quantized notation in the leads, because this 
notation readily connects to the formalism used later on. 



The stationary density matrix is found using the well- 
established and controlled Bloch-Redfield theoryii&. This 
is a systematic technique for deriving generalized mas- 
ter equations within Born approximation in Hi, eq. 
which includes all relevant non-markovian parts. This 
approach has been shownil to be numerically equivalent 
to formally exact path integral methods for the Spin- 
Boson model in the weak-coupling limit. The Redfield 
equationsi^ for the elements of the reduced density ma- 
trix p in the molecular basis read 
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where tOnm = {En — Em)/h are the appropriate energy 
splittings and Rnmki are the elements of the Redfield ten- 
sor. They are composed out of golden rule rates involving 
Hi from eq. n, m, k and I can be either -I- (molecu- 
lar excited state) or — (molecular ground state). The £'s 
are the eigenenergies of the two molecular states. Due 
to the lack of symmetry, this leads to a huge number of 
processes contributing to each termi^. We are only in- 
terested in stationary solutions here. A full treatment of 
the simple case with 7 = can be found in Rcf. 6. 
The current is derived from the standard formula— 



I(t) 
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where is the particle number operator on the left dot 
in the interaction representation and the transformed in- 
teraction Hamiltonian Hi from eq. JHl is also taken in the 
interaction picture. Carrying out the integration in equa- 
tion JSJ and rotating back to the Schrodinger picture, we 
get a time-independent expression for the current /. Us- 
ing the stationary occupation probability of the molecu- 
lar ground {p _st) or excited state (p++^st), we obtain 

for the expectation value of the stationary current 

1st = tr(psti") = p++,stl++ + P — ,st/ — , (9) 

where we find from balancing relaxation processes in the 
Bloch-Redfield equation, eq. (TJ 
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The current amplitudes I++ and / in eq. @ are 

of the same form as the contributions to the Redfield 
tensor. We emphasize, that the choice of processes from 
all possibilities is very distinct. As an example. Fig. [21 
displays a variety of possible processes in such a double- 
dot system. Processes of the type displayed in Fig. [3 

(a) contribute to the relaxation but do not carry current, 

(b) shows a process which carries current but does not 
relax the state, and (c) relaxes and carries current. The 
phase information of the quantum state is lost in all three 
pictures of Fig. |21 Consequently, one must not confuse 
cotunneling rates with relaxation rates. 
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FIG. 2: Examples for relevant processes: (a) a relaxation 
process without current, (b) current without relaxation (only 
dephasing) and (c) a process that carries current and also 
relaxes the system. 



We now turn to the discussion of the results. AU inter- 
nal energies £as and 7 are normahzed in units of the bias 
voltage V, the stationary current 1st in terms of Iq = eT. 
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FIG. 3: Stationary current 1st /lo for different Eaa/V as a func- 
tion of the coupling 7/I/. (with T = 140 mK, V = 5.170 /iV 
and /iav = (/ifl + /ii)/2 = 75.832 /jeV and F = 1 GHz.). 

In Fig.|3|the current at fixed bias voltage as a function 
of the inter-dot coupling is shown. The sign of Eas plays a 
role, as one can see above, for an intermediate Sas-regime. 
This effect is more pronounced in I{V), see Refii^. Close 
to 7 = 0, the curves all turn to zero, because at that 
point the dots are disconnected and no current can flow. 
However, a number of curves, the ones with £as/^ < 
1, exhibit an intermediate maximum at low 7 next to 
a very sharp minimum at 7 = 0, which sometimes is 
hardly resolved. At high j ^ V, the stationary current 
saturates into a value, which for our parameters turns 
out to be about /o,st/^o = 7.5 • 10~^. Remarkably, this 
is half the value of the current at the aforementioned 
low-7-Maximum. This is the central result of this paper. 

These regimes can be classified in terms of the level 
splitting 2S^^: At < 26, the energy V supplied from the 
leads is only sufficient to use one of the molecular states 
for transport (elastic cotunneling) whereas at V > 26, 
both states participate and also inelastic processes con- 
tribute, i.e. there is a second current channel, which car- 
ries the same contribution of Iq . The crossover naturally 

occurs at 7 = ~ ^"Lj which can only be reached if 

^as/V < 1/2. As long as 7 is not too low, the coupling 
to the leads is the limiting element for the current flow. 
Only if 7 < Eas, the double-dot eigenstates become local- 



ized and the inter-dot coupling becomes the current bot- 
tleneck. Consequently, associated dips have a half-width 
of Eas for low temperatures and bias voltages and can 
thus be extremely narrow. We would like to remark that 
the notion of transport "channels" is appropriate here, 
because cotunneling is a coherent transport process. 
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FIG. 4: Stationary current Ist/ht for different values of 'y/V as 
a function of the asymmetry energy Eus/V (with T — 140 mK, 
V = 5.170 pN and Mav = 75.832 /xeV and F = 1 GHz). 

Fig. 21 shows the dependence of the stationary current 
on E^jV . It confirms the interpretation of Fig. O The 
plot is only weakly asymmetric to e^^lV — 0. At zero 
asymmetry, £as/^ = 0; the condition for charge trans- 
port is ideal, \J ^ — maximum and therefore 

the current is only governed by the inter-dot coupling 
7/y, resulting in a zero-asymmetry maximum. 

Still, all three transport regimes can be recognized in 
Fig. 21 The ^ /V = curve shows, that the stationary 
current Ist/Io is exactly zero as expected. For growing, 
but small values of j/V, the maximum at Sa,s/V — 
reaches the highest value /st/^o — 2/o,st/-fo at about 1.5 • 
10~^ (like in Fig.EJ, corresponding to two open transport 
channels (elastic and inelastic). If we raise 7/V further, 
the height of the peak goes down again and saturates at 
Ist/Io = Io,st/Io ~ 7.5 • 10"'', corresponding to only the 
elastic channel being open. 

The three transport regimes are summarized in Fig.|Sl 
i) the atomic limit (no transport) 7 < Eas, ii) the 
two-channel case (inelastic cotunneling) Sas < 7 < 

and iii) the one-channel case (elastic cotun- 
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neling) j> ^^-el^. 

These conditions show that indeed cotunneling can be 
used as a tool for investigating the energy spectrum of an 
undisturbed artificial molecule'*. The cross-over between 
the elastic and the inelastic cotunneling in dependence of 
the applied bias voltage has recently been observedi^. A 
similar conclusion was found in Ref. I20I 

Although the notion of (elastic and inelastic) cotunnel- 
ing was already introduced very earljsii, its consequences 
for realistic quantum dot systems have only been dis- 
cussed very recentlySS,, along with detailed and accurate 
experiments on small semiconductor quantum dots^ii^ 
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FIG. 5: Limits for the three transport regimes with V = 
5.170 

becoming available. The sharp crossover between elas- 
tic and inelastic cotunneling, which we discuss, has been 
identified in a vertical quantum do^ by changing the 
transport voltage. Ref. i;^ studies cotunneling in a paral- 
lel double-dot topology, using again cotunneling and the 
elastic-to-inelastic crossover as a spectroscopic tool and 
tuning the inter-dot coupling in situ. In both cases, the 
narrow regime of decoupled dots would not have been 
accessible through a conductance measurement. Some 
of the experimental issues have been theoretically ad- 
dressed in Ref. I^^- In that case, however, the behavior 
of a single multilevel dot system was modelled with phe- 
nomenological couplings to the leads, whereas we take a 
realistic model and only by this manage to predict ef- 



fects which e.g. depend on the serial dot topology of the 
sample. Note, that parts of the double-dot literature fo- 
cus on phonon/photon assisted trans port (see e.g. Refs. 
l2lll2p for experiments and Refs. "23*2? for theory). Unlike 
Ref. I25L we concentrate on the Coulomb blockade regime 
and do not consider cotunneling at resonance. In Ref. 

a different approach to the problem was developed, 
in which the master equation is carried to second order 
instead of using a Schrieffer- Wolff transformation, and a 
few setups simpler than ours are studied. Our approach 
does not require the molecule to be artificial, in principle 
it can be applied to "real" molecules^. In contrast to 
the approach in Ref. it permits to take into account 
charging effects, however, the Schrieffer- Wolff transfor- 
mation is clearly a laborious step for larger systems. 

To conclude, we analyzed the stationary coherent co- 
tunneling current 1st through a double quantum dot sys- 
tem or artificial molecule. As a function of the inter- 
dot coupling strength it displays a rich, non-monotonic 
structure, which enables us to perform "molecular cotun- 
neling spectroscopy" . Strikingly, we have shown that at 
fixed bias voltage, the current is highest, if the dots are 
weakly-to-intermediately connected, such that the inter- 
dot coupling is at least as strong as the coupling to the 
leads, but the splitting of the molecular wavcfunctions is 
still smaller than the transport voltage. 
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